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of narrow band frequencies across a wide spectrum 
are used in a frequency hopping transmission pro- 
tocol, the accuracy in determining the phase differ- 
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DIFFERENTIAL RANGING FOR A FREQUENCY-HOPPED 
REMOTE POSITION DETERMINATION SYSTEM 

FIELD OF THE INVENTION 
5 The present patent application relates to radio communication 

and transmitter location system. 

BACKGROUND OF THE INVENTION 
Location of radio transmitters is a technology that has been 

10 used for many years and in many applications. For radio location, the use of 
mobile transmitters and a plurality of Fixed base station receivers is generally 
known in the art where the fixed receivers locate the source of the radio signal 
by triangulation of the radio signals. Triangulation is a technique where a 
plurality of fixed base station receivers each determine the approximate 

15 direction of the radio signal from the transmitter. The approximate direction 
lines are then drawn on a map and the location of the transmitter is defined at 
the cross point of these lines. Since the direction lines arc approximate, and 
typically three fixed base station receivers are used, the intersecting area 
between the three lines is usually shaped like a triangle since the lines rarely 

20 meet at a single point. The transmitter is then located within the triangle. 

Triangulation of radio signals is a costly undertaking since at 
least two receiving antenna stations must be situated at different diverse 
geographic locations to perform even rudimentary triangulation. Effective 
triangulation is more realistically achieved with three receiving locations and 

25 accuracy is increased even more with four receiving towers. The more 
receiving stations constructed, the better the accuracy in locating the 
transmitter, but also the great the cost. 

One such direction finding technique is described in US Patent 
No. 5.379,047 entitled "REMOTE POSITION DETERMINATION SYSTEM" 

30 issued January 3, 1995. In that patent, the direction of the incoming radio 
signal is determined by using an array of receiving antennas to measure the 
phase difference between the antennas of the array and thereby calculate the 
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incident direction of the radio carrier signal. Although this patent has many 
novel and useful advantages, it requires that each base station use an array of 
antennas and at least two and typically three bases. Quite often the cost of 
building and maintaining these base stations can be quite high, especially if 
5 they are located in urban areas where roof-top space must be leased and the 
cost of the lease is determined by the number of antennas placed thereon. 

Another disadvantage of triangulalion for location of radio 
transmitters is the construction of such a system in sparsely-populated areas. 
In the example of a highway across a desert area, it is cost prohibitive to 

10 install a large number of base stations to track a narrow stretch of road. The 
base stations would have to be placed on both sides of the road and spaced 
quite far back from the road. 

There is a need in the art to locate a radio transmitter using at 
least three and in special cases four one-antenna receivers. A receiver with 

15 one antenna will minimize the construction and lease costs of a location 
system. There is a further need in the art to determine the location of low 
power, frequency hopped spread spectrum radio transmitter. 

SUMMARY OF THE INVENTION 
20 The present invention solves the above-mentioned problems 

with locating transmitters and other problems which will be understood by 
those skilled in the art upon reading and understanding the present 
specification. The present invention uses a modified time-of-arrival technique 
to determine the location cjf a frequency hopped spread spectrum radio signal. 
25 The transmitter simultaneously transmits two radio frequency carriers having 
different frequencies such that a phase difference is observed between the two 
carriers at a distance from the transmitter. The phase difference is 
proportional to the range from the transmitter that the carrier signals are 
observed. The two carrier signals from the single transmitter are received by 
30 at least three and in special cases four base stations which calculate the 
differential time of arrival based on the phase differences of the received 
carriers. The calculated phase differences are then sent to a central location 
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which locales the position oflho transmitter based upon a planer hyperbolic 
location algorithm. Since a large number of narrow band frequencies across a 
wide spectrum are used in a frequency hopping transmission protocol, the 
accuracy in determining the phase differences between the two carriers is 
5 increased and the immunity to interference is also increased. 

BRIEF DESCRIPTION OF THE DRAWINGS 
In the drawings, where like numerals refer to like components 
throughout the several views. 
10 Figure 1 is a diagram of an integrated frequency hopping spread 

spectrum communication and location system; 

Figure 2 is a block diagram of a reverse paging and mobile 
location system using the infrastructure of a standard paging system; 

Figure 3 describes the synchronization and message format of 
1 5 the outgoing paging signals from the base stations; 

Figure 4 describes the format of the frequency-hopped spread 
spectrum signal transmitted by the remote mobile units; 

Figure 5 is a hyperbolic plot of the differential ranging of the 
present location system; 
20 Figure 6 is a diagram of the differential time-of-arrival aspects 

of the present location system; 

Figure 7 shows the computations used to perform location using 
differential time-of-arrival ranging; 

Figure 8 describes in more detail the principle of transmitting 
25 two simultaneous frequencies to two receiving base stations; 

Figure 9 is block diagram of the differential ranging transmitter: 
Figure 10 is a diagram of the IF circuit of the differential 
ranging receiver in the base station which uses two separate convenors: and 

Figure 1 1 is a phase plot showing the effect of sources of errors 
30 on the differential ranging technique. 
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DKTAILHI) DI-SCRIPTION OK Till-: PRm'KRRlH) KMHODIMHNT 
In the following detailed description of the preferred 
embodiment, reference is made to the accompanying drawings which form a 
part hereof, and in which is shown by way of illustration specific preferred 
5 embodiments in which the inventions may be practiced. These embodiments 
are described in sufficient detail to enable those skilled in the art to practice 
the invention, and it is to be understood that other embodiments may be 
utilized and that structural, logical and electrical changes may be made 
without departing from the scope of the present inventions. The 

10 following detailed description is, therefore, not to be taken in a limiting sense, 
and the scope of the present inventions is defined by the appended 
claims. 

Svstem Overview 
In the preferred embodiment of the present invention, the base 

15 stations and central site are similar to the base stations 

and central site described in U.S. Patent Number 5,430,759 entitled "LOW- 
POWER FREQUENCY-HOPPED SPREAD SPECTRUM REVERSE 
PAGING SYSTEM" issued July 4, 1995. and the processing structure of this 
system operates similarly to the present system. The present invention allows 

20 the accurate location of the frequency hopped, spread spectrum transmitters 
even in noisy and multi-path environments such as those found in urban 
environments and irregular terrain environments. 

The present invention is directed to a location system for use in 
locating remote mobile units (RMU's). The remote mobile units may be 

25 located in motor vehicles, located on the person of people, carried within 
containers or packages, or any number or variety of mobile carriers. The 
location system comprises a plurality of base stations (at least three and in 
special cases four) which can transmit signals to the remote mobile units and 
receive signals from the remote mobile units. The base stations may send 

30 messages and information to the remote mobile units, the base stations may 
send alert or interrogation commands to activate the remote mobile units, or 
the remote mobile units may be activated locally without any prompting. 
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Thus the remote mobile units may be activated locally (lor example by the 
wearer or vehicle) or remotely by the base station. The remote mobile units 
may be used by a person as a reverse pager, as an emergency locator or as a 
communication device. The remote mobile units may also be used as part of 
5 a vehicle to locate the vehicle if stolen. Tor vehicle tracking or as an integral 
communication device. The remote mobile unit may also be activated in the 
case of an accident by a deceleration switch. ABS system, air-bag or other 
crash-activated switch to automatically summon help. 

The remote mobile units receive messages as a standard paging 

10 device over licensed airwaves using a standard paging infrastructure. The 
remote mobile units, when activated, transmit low-power (less than one watt), 
frequency-hopped, spread-spectrum communication signals. The transmitted 
signals from the remote mobile units are received by the base stations which 
are then used to locate the remote mobile units. The base stations are 

15 equipped with a single receiving antenna. The location of the remote mobile 
units is determined by differential ranging using phase difference comparisons 
between two carrier signals simultaneously transmitted by the remote mobile 
unit. Unique algorithms are performed on the received signals at the base 
stations to eliminate the ambiguity caused by signal interference, phase 

20 ambiguity caused by errors and noise and multi-path reflections. 

The frequency hopping transmitter of the remote mobile unit 
produces two sine waves for each frequency hop. The phase difference 
between the two sine waves will vary with distance from the transmitter which 
is proportional to the range from the transmitter. There is an ambiguity factor 

25 in the range determination due to the phase differences between the two sine 
waves. The ambiguity is minimized by transmitting at least 53 hops per 
transmitted message so that there is a different phase difference within each 
frequency hop. Each hop change produces a step in frequency such that 53 
steps of various step distances in frequency are produced. 

30 At the receiver, the frequency hopping carriers are 

downconverted. sampled and digitized for processing in the digital domain. A 
range calculation is determined by comparing the phase differences of the two 
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sine waves of the two frequencies within each hop. The ambiguities are 
minimized by comparing the phase differences across the 53 hops in 
synchrony with the hopping protocol. The differential range between two 
bases is determined by collecting the phase differences from all the hops in 
5 two bases and performing in the central station an Inverse Fourier Transform 
(I FT) which will produce a single value for the differential range between the 
two bases to the transmitter. 

Each base station is equipped with a very accurate rubidium 
atomic clock so that the exact time of arrival of the frequency hopping signals 

10 is determined with reference to a very accurate standard time. Those skilled 
in the art will readily recognize that any highly-accurate clock means may be 
substituted for the preferred rubidium clock such as a GPS receiver. The 
central station compares the time of arrivals with t.,e knowledge that each 
base station has made its calculations based on the same time standard. The 

15 time of arrival information is then used to plot possible locations of the 
transmitter using a hyperbolic plot. The comparison of the time of arrivals 
from two base stations will produce a hyperbolic plot where the base stations 
are the foci. By plotting the time of arrival numbers on the hyperbolic plot 
and crossing the hyperbolic plots between two or more base stations, the exact 

20 location of the transmitter can be determined. 

To assist in calibrating the differential ranging system with the 
local base stations, a transmitter having a known location is used to 
periodically transmit a signal which is then used to synchronize the base 
stations. The calibration signals is used to adjust the local atomic clocks of 

25 the base stations. It is possible to use less accurate clocks in the present 
system by increasing the number of calibrations performed in a day. In a 
calibrated system the accuracy of location has been measured at 26 meters at 
minimum received signal strength of -138dbm when operating within the 
parameters described below. 

30 The problems associated with multipath reflections and fading 

are known in the art. The multi-path problems and fading problems are 
usually frequency dependent so the use of a large number of frequency hops 
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virtually eliminates these problems. The ambiguity problems are also solved 
in a fashion similar lo thai described in the aforementioned copending patent 
application. 

5 Summary of I Ivperbolic Location with Differential Time of Arrival 

The operation of the present differential ranging system may 
also be termed differential time-of-arrival (TOA) since the differences in the 
arrival limes at ihc base stations are used to determine the position of the 
receiver. As is described in more detail below, the TOA of the two carrier 

10 frequency signals for each remote mobile unit are estimated at each base 

station for each hop. The TOA is proportional to the range R of the signals. 

The operation of the present system is basically performed in 
two major steps. In the first step, estimated range R„ is determined for each 
base station BS n , where n is the base station number. Due to the inaccuracy 

15 of the clock in the transmitter, the range by itself is not an accurate indicator 
of the distance of the remote mobile unit from the transmitter. Since the base 
stations are equipped with very accurate clocks which are periodically 
synchronized to a common time base, the base stations can report the exact 
time of arrival of the signals. 

20 In the second major step, the lime of arrival information from 

each base station and the range information (calculated across the multiple 
hops using an Inverse Fourier Transform) are sent to a central station where a 
range differential between two base stations is plotted, as shown in Figure 5. 
The range differential forms a hyperbolic plot, for example 501a and 501b, 

25 between base station BS, and base station BS 3 . We know that the range R, is 
greater than the range R 2 so that the hyperbolic plot 501b in dashed lines in 
Figure 5 can be ignored. A second hyperbolic plot 502a and 502b is formed 
between base station BS 2 and base station BS 3 . The intersection point 
between the two hyperbolic plots is the location of the remote mobile unit 

30 RMU. 
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Comparison of Differential Raiminti to (jPS 
The operation of Ihc present differential ranging system is quite 
different from the GPS (Cilohal Positioning Satellite) system. The UPS is a 
navigation system that provides location for the vehicle that carries it. In 
5 order to monitor the vehicles, the system requires a communication system for 
delivering the location to a central office. The navigation process of the GPS 
is virtually the inverse process of the dua! frequency differential ranging 
system of the present invention. The theoretical accuracy of the present 
invention is only dependent on the bandwidth, signal to noise ratio and some 
10 ambiguities described below. Although commercial GPS systems are made 
inaccurate on purpose, a comparison can still be made. 

There are several facts in favor of the present differential 
ranging system for the task of remote mobile unit location over GPS. The 
GPS system requires direct line of sight to the satellites, since it has a very 
15 sensitive receiver. Differential ranging is transmitting in the 900 MHz band 
and thus does not require line of sight: it could operate from within buildings. 
Differential ranging does not require any other device other than the low cost 
RMU on the vehicle which makes it significantly lower cost. The power 
consumption of the differential ranging system is significantly lower since it 
20 requires only a low power transceiver, thus it is operates on batteries for a 
long time (over 1 month on commercial batteries). The differential ranging 
system operates using frequency hopping and is less susceptible to multipath 
and interference than the GPS which operates in direct sequence spread 
spectrum. From measurements performed on the present direction finding 
25 system, it is by far less susceptible to ground multipath than the GPS. The 
complex receivers of the system are only in the base stations. The 

complex receivers of the GPS are in every vehicle. 

Frequency Hopping System Overview 
30 In the preferred embodiment of the present invention, 

transmitters and base stations similar to those used in the present invention are 
described in the aforementioned copending U.S. Patent Number 5.430,759 
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entitled "LOW-POWKR l-'RI-OlJIiNCY-IIOPIM-D SPRHAI) SPKCTUUM 
RKVKKSK PAGING SYSTKM" issued July 4. IW5. The infrastructure of 
this reverse paging system serves as the hasis of the present and the 

reverse pagers of this system operate similarly in the present system using a 
5 single carrier frequency instead of the two carrier signals of the present 
invention. The system allows the accurate location of the reverse 

pagers even in noisy and multi-path environments such as those found in 
urban environments and irregular terrain environments. 

The environment in which the present system is operated is 

10 shown in Figure 1. The system is designed to operate in a 

communications environment known as NexNet™ in which a plurality of base 
stations service a particular area using the NexNet™ communications system. 
Each remote mobile unit receives standard paging messages as a downlink and 
can provide talk back messages as an uplink channel using frequency hopping 

1 5 spread spectrum radio frequency communication. The paging transmitter of 
the base stations 1 13 typically transmits with 200 Watts of power to send 
messages using a standard paging formal such as POCSAG or other industry 
standard formats. Each remote mobile unit contains a paging receiver with a 
sensitivity of -115 dBm. a frequency-hopping dual frequency transmitter with 

20 a transmit power of one Watt and double side band modulator. Each base 
station includes a receiving antenna for the frequency hopped carrier signals 
and special processors for determining the phase difference between the two 
carriers of the double side band modulation. 

Figure 2 depicts the major components of the two-way paging 

25 system in the aforementioned U.S. Patent Number 5,430,759 entitled "LOW- 
POWER FREQUENCY-HOPPED SPREAD SPECTRUM REVERSE 
PAGING SYSTEM" issued July 4, 1995. The reverse paging terminal 1 10 at 
the central site operates to provide synchronization and messaging information 
through the paging terminal 1 16 to the reverse pagers 100 (also known as 

30 remote mobile units 100) via direct links to the base stations BS,, BS 2 and 
BS^ through ground based radio links (not shown) or through a satellite 
uplink/downlink using a geostationary satellite (not shown). The base stations 
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BS,. BS ? and BS, include transmit and receive towers 1 13a. 1 13b and 1 13c. 
respectively and base station terminals 201a. 201b and 201 c. respectively. 
Terminals 20la-201c arc required for producing the accurate synchronization . 
information needed to be transmitted to the remote mobile units and for local 
5 processing of the received messages for transmitter location. This 

synchronization information is used to coordinate the frequency hopping 
transmissions and to coordinate the response of messaging from the plurality 
remote mobile units 100 so as to minimize collisions within groups of remote 
mobile units and eliminate collisions between groups of remote mobile units. 

10 

Base to Remote Synchronization 
Standard paging messages sent from the base stations BS,. BS, 
and BS*, to the plurality of remote mobile units 100 are, in the preferred 
embodiment, sent as digital data encoded in the POCSAG paging standard. 

15 These messages may be used to interrogate the remote mobile units to activate 
the remote mobile unit to allow the base stations to begin the location process. 
Typically the paging channel has a center frequency of 143.160 MHz, with an 
NRZ FSK data rate of 512 bps or 1200 bps. Other bit rates such as 2400 
baud (bps) are also feasible. Figure 3 describes the POCSAG paging 

20 communications protocol as modified for use by the preferred embodiments of 
the present invention. In the top line of Figure 3, a greatly compressed time 
line of digital data transmitted according to the POCSAG protocol is shown. 
Batches of messages are transmitted in groups as shown in the details in the 
subsequent lines below the top line of Figure 3. In the second line of Figure 

25 3. a 1.0625 second interval (for 512 baud) is shown in which 544 bits are 
transmitted as a single batch. The batch is preceded by a synchronization 
code word SC as shown in the third line of Figure 3. 

The synchronization code word within each batch is followed 
by eight frames of digital data. Each frame is divided into two portions, an 

30 address portion and a message portion. The address code word of the 

message of frame 2 of Figure 3 is shown in line 4 while the message code 
word of the second half of frame 2 is shown in line 5. The address code 
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word is preceded by a digital zero followed by IK address bits, two function 
bits and 10 check bits. The address code word is followed by an even parity 
bit. The message code word portion of the frame is preceded by u digital one 
followed by 20 message bits which are followed by 10 check bits and a single 
5 even parity bit. Thus each frame is comprised of 64 bits divided into two 32 
bit sections. 

Synchronization of the base station terminal 200 and the remote 
mobile units 100 is necessary to ensure the units 100 are transmitting at the 
same time that the base stations are listening. Synchronization is also 

10 necessary to coordinate the division of the large number of remote mobile 
units into groups so that members of one group use different frequency 
hopping patterns from members of other groups. Synchronization of the 
remote mobile units 100 is accomplished by inserting a special frame into the 
POCSAG data which is used to synchronize the units. 

1 5 The purpose of synchronization between the reverse paging 

terminal 1 10 and the remote mobile units 100 is to determine where along the 
pseudo random noise code the frequency hops are to be followed and to 
determine the exact times for transmitting frequencies from within any of the 
hops. This also enables the dynamic changing of a remote mobile units group 

20 membership such that if one group is experiencing a large number of 

collisions due to simultaneous transmissions, the reverse paging terminal 1 10 
may re-allocate some of the remote mobile units within that group to new 
groups to minimize collisions. 

Referring once again to Figure 3. eight frames of information 

25 are transmitted in each burst using the POCSAG format. Remote mobile units 
100 may be assigned to a specific frame within the transmission so that the 
remote mobile units, once recognizing the synchronization code word, can 
scan a specific frame for that remote mobile unit's address. Once the address 
is found, the remote mobile unit can determine any group changes that may be 

30 required to re-allocate that remote mobile unit to a different group. In 

addition, the POCSAG format is used to transmit a fine time synchronization 
code. The fine synchronization code is a transmission of a time pulse at an 
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exact time synchronized to a GPS (Global Positioning System) clock to 
synchronize all the remote mobile units 100 Tor time of transmission, l or 
example, periodically during the day the reverse paging terminal will send a 
synchronization code within the POCSAG code word which is sent at a very 
5 precise time. In order to ensure that a precise time pulse is sent, the reverse 
paging terminal 1 10 receives accurate time information using a GPS antenna 
to receive accurate time of day information. The time used to send the 
synchronization pulse is when the day clock reaches exactly some multiple of 
0.9 seconds in the preferred embodiment. In this synchronization information. 
10 20 bits of information are transmitted to give the accurate time of day 
information. 

In each of the remote mobile units 100. the microprocessor 
compares this accurate time pulse which will indicate the exact time of day 
and compare it to its own day clock. The clock within each microprocessor is 

15 accurate down to a few milliseconds, but the lime at which the 

synchronization pulse occurs should have a resolution much finer than that 
such as down to 0.1 milliseconds for time of day. In this fashion, each of the 
microprocessors in each of the remote reverse paging devices can periodically 
realign its day clock to know within a millisecond the exact time. Each 

20 microprocessor does not actually realign its clock but changes a clock offset 
within memory so that it understands how far off its own internal clock is and 
can make the adjustment when using that clock to determine when to start 
transmitting information. 

The synchronization pulse is only transmitted every few 

25 minutes. However, the resolution of the start of the message indicating the 
synchronization pulse is very accurate, it being transmitted at 0.090000 
seconds GPS time after a fixed time of day. such as 12:00 GMT. This GPS 
time is accurate to at least within 100 nanoseconds. 

An overview of the transmission formal of the remote mobile 

30 unit is shown in Figure 4. The actual transmission of information from the 
remote mobile units 100 is done using Differential Bi-Phase Shift Keying 
(DBPSK) modulation on a frequency hopped carrier of less than one watt. 
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The transmission of information from the remote mobile units 100 on the 
frequency hopped carrier may also be done using Frequency Shift Keying 
(FSK) modulation. The same information is modulated on the two carriers in 
the preferred embodiment of the present invention allowing the base station to 
5 decode the information from either carrier or both carriers. The dual carriers 
arc preferably two different narrow hop frequencies. 

Typically a single transmission consists of 53 hops or 53 
changed frequencies selected from a list of narrow band frequencies. The 
frequency selection is based on a pseudo-random noise code list pointing to 

10 the frequency selection list. The synchronization information tells the remote 
mobile unit 100 where along the pseudo random noise code it should be 
synchronized for transmission of its message and tells exactly the time of day 
so that the remote mobile unit 100 knows exactly when to start transmitting 
the specific frequency so that the base stations BS1. BS2 and BS3. labeled 

15 1 13a, 1 13b and 1 13c, respectively, are looking for that frequency at the same 
time. 

In operation. 200 frequencies are used by the remote mobile 
unit 100 and the base stations and internally stored in a list numbered Fl 
through F200. For a specific message. 53 frequencies will be used to transmit 

20 the entire message. These 53 frequencies are selected based on a 1.000 
member pseudo-random noise code. 

The use of the accurate synchronization signal periodically 
broadcast via the outbound paging signal enables the remote mobile units to 
use lower accuracy components thus reducing the manufacturing cost of 

25 remote mobile units. For example, high accuracy crystals to track the time of 
day within the microprocessor are available with an accuracy of three parts 
per million. Thus, a time drift of approximately three micro seconds per 
second or 180 microseconds in a minute is the known drift. There are also 
time inaccuracies which are introduced due to variable path length between 

30 the paging tower to the remote mobile unit. By employing crystals which are 
cheaper and have an accuracy of the order 50 parts per million, the amount of 
time-of-day drift normally wouldn't be tolerable. However, by using the 
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synchronization information transmitted on a regular basis from the reverse 
paging terminal, the microprocessor can continually correct its own internal 
day clock so that accurate time of day measurements are always maintained. 
The microprocessor estimates the momentary inaccuracy of the crystal by 
5 tracking the drift across several synchronization transmissions and dynamical ly 
adjusts for the frequency drift of the crystal and the offset using internal offset 
registers for accurate time of day information. 

Counters are employed within each microprocessor of the 
reverse paging units to compensate for the offset of the frequency based on 

10 the synchronization time information. There are generally two major factors 
which affect the drift in a crystal: temperature and acceleration. Most of the 
drift is due to temperature, and the remaining drift components are negligible. 
The frequency drift in a crystal due to temperature is very slow, on the order 
of 50 Hz over 10 seconds. During a single day the temperature can change by 

15 20 or 30 degrees fahrenheit, requiring a time update from the GPS clock 
approximately every five minutes. 

Remote Mobile Unit Transmission Format 
The signal sent from the remote mobile unit 100 to the base 
20 stations is a spread-spectrum, frequency-hopped transmission using differential 
bi-phase shift keying (DBPSK) modulation on the frequency-hopped carrier to 
transmit digital information. The transmission of information from the remote 
mobile units 100 on the frequency hopped carrier may also be done using 
Frequency Shift Keying (FSK) modulation. The frequency hops are relatively 
25 slow, the frequencies transmitted are very narrow and the transmission power 
is extremely small. 

The maximum peak output power of transmission from remote 
mobile unit 100 is limited to less than one Watt to allow use of the 902-928 
MHz ISM band in the United States without the need for licensing the remote 
30 mobile units as allowed by FCC regulations defined in 47 C.F.R. §15.247. 
Those skilled in the art will readily recognize that other frequency bands and 
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transmissions power levels may be employed depending upon l-VC licensing 
requirements or oilier frequency licensing requirements of other nationalities. 

The use of an accurate crystal u> control each frequency of 
transmission is required within each remote mobile unit 100. l ; or example. 
5 high accuracy crystals to transmit the narrow bandwidth frequencies used for 
the frequency hopped transmissions are available with an accuracy of three 
parts per million. At 900 MHz. a 3 ppm drift would place a single frequency 
somewhere within a 2.7 KHz band. To tolerate frequency drift due to aging 
and temperature, each individual frequency of the frequency hopped signal is 

10 allocated to a 7.5 KHz band or channel, even though the actual frequency is 
on the order of 200 Hz wide skirt within this 7.5 KHz allocated bandwidth. 
Those skilled in the art will readily recognize that by using alternate 
components, the frequency channels (individual frequency of the frequency 
hopped signals) of 7.5 KHz allocated bandwidth may be wider or more 

15 narrow depending upon the overall allocated bandwidth for the system. For 
example, 1 KHz or less bands may alternatively be allocated per channel. 

Tests on this system have shown that by processing the 
received signals at the base stations entirely in the digital domain using the 
combination of unique Fast Fourier Transform algorithms 

20 to locate and retrieve the frequency hops and by using a 

combination of unique confidence algorithms with a plurality of error 
correction codes, the receiving base station is able to pull the response 
information from a very low power signal from a distance of up to 45 
kilometers (28 miles) in a flat terrain. In a rather noisy urban environment, a 

25 range of 24 kilometers (15 miles) is the norm. The information within the 
signals is accurately decoded even in severe multipath and noise conditions. 

As shown in Table I, the remote mobile unit message format 
consists of a preamble and the message body spanning a total of 53 frequency 
hops. Those skilled in the art will readily recognize that longer messages may 

30 be transmitted using the preferred embodiment of the present invention, and 
the messages format described here is illustrative but not limiting. Much 
longer message hops lo transfer more digital data is also implemented but not 
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described here. Of course, Ihosc .skilled in ihe ari will readily recognize thai 
shorter messages than those described below are equally possible lor the 
preferred embodiments of ihc present invention. The message length and 
number of transmission bops arc a mailer of design choice. 
5 The message preamble consists of a predefined code of ones 

and zeros lo gel Ihc attention of the base unit receiver to pull the message out 
of the noise. The preamble consists of 165 bits transmitted across 5 hops, that 
is. transmitted using DBPSK (Differential Bi-Phasc Shift Keying) or 
Frequency Shift Keying (FSK) on five different frequencies selected from the 
10 frequency list with the specific frequencies selected based on the PN (Pseudo- 
random Noise) Code list stored within the remote mobile unit. The sequence 
location within the PN code that the remote mobile unit will begin to follow 
is based on the synchronized time of day. Within a single bop (a single 
carrier frequency), the carrier phase is modulated 33 times to encode the 
1 5 predefined one-zero pattern of the preamble. 

The message body follows the preamble and consists of three 
groups of data. Each group consists of 30 actual data bits so that the entire 
message is, in the preferred embodiment of the present invention, 90 total data 
bits (although other bit length messages may be chosen). The actual data 
20 encoded within these 90 bits is described above and may be in any convenient 
coded format. Those skilled in the art will readily recognize that a wide 
variety of message formats and encoding of the data bits may be used without 
departing from the scope of the present invention. The encoding 

described here, however, has been proven effective in retrieving the data bits 
25 buried in background noise with a high degree of accuracy and a low actual 
error rate. 

TABLE 1: Remote Mobile Unit Message Format 
Preamble is 165 bits (33 bits x 5 hops) 
Message is 48*33 transmitted bits 
30 (Message is 90 bits actual data) 

[] = One Frequency Hop 



35 



T 1 1 1 1 1 1 1 1 1 I I I I I 1 1 1 I J. 1 r 



preamble J Message Body Spread Over 4 8 Hops 
, , i i | i i i i i » i » i i I I i i — I — i — l— - — J — 1 
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Outer Messauc Codinu 
r-ach of the three groups of message daia (30 bits each) arc 
1 KM I encoded using a standard 30.63 BCII code and with a single parity bit 
added to form a 64-bit word. This encoding decreases the error rate from 10° 
5 to 10' 5 . This encoding, documented and understood by those skilled in the an. 
can correct up to 6 errors or delect up to 13 errors. Detection of corruption of 
a data word that cannot be reconstructed will cause the base to request a 
second transmission of the acknowledgement message. 

10 Inner Codiniz and Interleaving 

The inner coding of the message will protect the integrity of the 
message with an error rate as high as 25%. Each block of 64 bits of data 
(corresponding to a groups of 30 bits and earlier encoded by a standard 30.63 
BCH code) is split into two sub-blocks of 32 bits (sub-blocks A and B of 

15 Table 2). and a reference bit is added to each sub-block to assist the 
differential encoding to provide a reference bit to the DBPSK or FSK 
decoder. The 33 bit sub-blocks are transmitted over one frequency hop each 
and are replicated 8 times so that the 64-bit block traverses 16 frequency 
hops. In transmission, the 33 bit sub-blocks are interleaved to further reduce 

20 loss of data, as shown in Table 3, where sub-blocks A and B of Table 2 

correspond to the first group of 30 bits, sub-blocks C and D. correspond to the 
second group of 30 bits. etc. The total message is 53 hops where each hop is 
180 msec in length making the duration of a single message 9.54 seconds. 
TABLE 2: Interleaving Format for Sub-block 

25 

A = 1 reference bit and 32 data bits = 33 bits 
B = 1 reference bit and 32 data bits = 33 bits 

[] = One Frequency Hop 

30 



A 


B 


A 


B 


A 


B 


A 


B 


A 


B 


A 


B 


A 


B 


A 


B 
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TABLE 3: Inner Coding and Interleaving of Sub-blocks 

A = first 33 bits of lsz block 

B = second 33 bits of Isl block 

C = first 33 bits of 2nd block 

D = second 33 bits of 2nd block 

E = first 33 bits of 3rd block 

F = second 33 bits of 3rd block 

[] = One Frequency Hop 



15 



20 



— Till 

preamble 
(ill. 


A 


B 


A. . .B 


A 


B 


C 


D 


C ... D 
1 


C 


D 


E 


f| . 


E 


F 


E 


F 



Those skilled in the an will readily recognize that a wide 
variety of data interleaving may be utilized to effect beuer error tolerance and 
may be substituted for the interleaving described here. Such alternate 
substitute interleaving means are CIRC (Cross Interleaved Reed Solomon 
Code) used in CD (Compact Disc) recording media operating either at the 
block level or at the bit level. 



Sintile Hon Format 
The acknowledgment signals are transmitted by the remote 

25 mobile units 100 in a 1.5 MHz band selected from within the 902-928 MHz 
spectrum. The 1.5 Mhz band is divided into 7.5 KHz channels to provide 200 
channels available in which the frequency hops can occur. Thus, each 
frequency hop is a channel 7.5 KHz wide in which a carrier frequency is 
transmitted. For example, channel one will have a frequency Fl at 

30 902.00375 MHz + -3.75 KHz, channel two will have its center carrier 
frequency at 902.01025 MHz + -3.75 KHz. etc. 

Each transmit frequency of each hop will thus be centered at 
the approximate mid-point of the assigned channel band: however, due to 
inaccuracies in the remote mobile unit circuits and reference crystals, the 

35 actual transmit frequencies will vary between units. If high quality crystals 
are used to accurately produce the required frequencies, very little drift off the 
center frequency will result. In the preferred embodiment of the present 
invention, low cost crystals are purposely employed to keep the per-unit 
manufacturing costs down. This will allow for a lower-cost product sold to 
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the user which will increase market penetration. Thus, reference crystals are 
preferred which have a frequency accuracy of 3 ppm such thai at 900 MHz. 
the statistical drift would be approximately 2700 Hz. The crystals center 
frequency within its nominal accuracy also drills due to aging and temperature 

5 variations, but this drift is slow compared to the transmission limes so the 
drift during a single transmission due to these latter variants is unimportant. 

A single frequency hop is shown in Table 4. The 15 
millisecond guard time preceding each hop is primarily a settling lime for the 
oscillator circuits of the remote mobile units to allow the internal oscillator 

10 circuit to lock onto the new frequency between hops. Each hop is transmitted 
at a single frequency in which the phase of the carrier is either at 0 degrees 
phase or 180 degrees phase in reference to the phase of the reference bit 
immediately following the quiet or guard time. Thus the first bit is a phase 
reference bit followed by 32 data bits exhibiting either zero phase shift or 180 

15 degree phase shift to encode the data bits as DBPSK (Differential Bi-Phase 
Shift Keying). In an alternative implementation, each frequency hop may be 
modulated using Frequency Shift Keying (FSK) in which two frequencies are 
used to transmit data bits. One hop frequency may indicate a logical one 
while a second hop frequency may indicate a logical zero. The frequency 

20 shift is minor and the frequency differential is contained within a single hop 
channel. 

Each bit of DBPSK or FSK is a transmission of 5 milliseconds 
of the hop carrier frequency either in phase with the reference bit transmission 
or 1 80 degrees out of phase. 
25 TABLE 4 : Single Frequency Hop Format 

Guard time (quiet) = 15 ms 

Single Bit = 5 ms of carrier DBPSK / FSK 

33 Bits plus guard time = 180 ms 

30 
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5 ms 
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5ms 
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2nd 


3rd 
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32nd 


35 


Time . 
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Bit 


Bit 
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Frequency Honninu Sequence 
All oflhc remote mobile units in the market serviced by the reverse 
paying terminal for message or location finding use the same pseudo random noise 
code to determine the frequency hops. The pseudo random noise code is a digital 
5 code of 1000 unique numbers. In the preferred embodiment of the present 

invention, the pseudo random noise code is stored in memory of each of the remote 
mobile units. Those skilled in the art will readily recognize, however, that a linear 
feedback shift register could be used to generate the pseudo random noise code on 
a real-time basis instead of using a look-up table which is presently in the preferred 
10 embodiment. 

The PN (pseudo-random noise) code list is stored in memory and 
maps to a frequency list. In the preferred embodiment of the present invention, the 
PN code list has 1,000 entries which repeat as a sequence. The control means of 
the reverse paging units continuously maintain a count of the proper location within 

15 this list. As described below, the time of day for all remote mobile units in the 
market served by the base terminal are periodically synchronized to ensure 
acknowledgment messages are synchronized to transmit the hop frequency at the 
proper time and to synchronize the location within the PN code list that each 
remote mobile unit will use to transmit. 

20 The 1,000 member PN code list maps to a 200 member frequency 

list. In order to allow a large number of remote mobile units to simultaneously 
operate in the same geographic market, the remote mobile units are divided into 
groups and the groups are assigned different sequence segment locations in the 
same 1.000 member PN list. Thus a remote mobile unit from group one will begin 

25 transmitting a hop at a frequency determined from a first location with the PN 
code, while a remote mobile unit from group two may begin transmitting a hop at 
a frequency determined from a second location in the PN code. The remote mobile 
units from group one and group two will complete their respective 
acknowledgement messages in 53 hops. Preferably, the sequence of the PN code 

30 used to determine the frequencies of the 53 hops for the remote mobile unit of the 
first group will not overlap the sequence of the PN code used to determine the 
frequencies of the 53 hops for the remote mobile unit of the second group. N4ore 
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preferable, the frequencies chosen based on ihe non-overlapping segments of ihe 
PN code list are orthogonal such that the same frequency is never used by two 
remote mobile units belonging to different groups. 

In the preferred implementation, the 1.000 member PN code list is 

5 divided into 160 hopping sequences. The remote paging units are divided into 40 
uroups with the members of each group synchronized to track the same location in 
the PN code list. The microcontroller of each remote mobile unit, regardless of its 
group membership, continuously runs through the repeating PN code sequence to 
stav in synchronization with the base unit and all other remote mobile units. Each 

10 group of remote mobile units is further divided into four subgroups such that the 
remote mobile units within each subgroup are assigned one sequence within the PN 
code list. Although the 53 hop sequence needed for each acknowledgement 
transmission may overlap the 53 hop sequence used by a remote mobile unit in 
another subgroup, the transmission sequences of a remote mobile unit of one group 

15 is chosen to not overlap the 53 hop sequence used by a remote mobile unit in 
another group. 

Base Station Design 
The analysis of the dual carriers and the decoding of the signals 

20 received by the base stations from the remote mobile units is done almost entirely 
in the digital domain. The carrier frequencies of the frequency hops are down- 
converted to a lower frequency in each base station and are then digitally sampled. 
The digital samples are then processed to locate the phase information of interest 
for direction finding and message decoding. Each base station is constructed with 

25 a plurality of digital signal processing based receivers which enable simultaneous 
message decoding and direction finding of a plurality of simultaneously 
transmitting remote mobile units. The dual carriers arc analyzed in parallel to 
determine the exact phase separation. Each bases station is also equipped with a 
rubidium atomic clock for accurate time-of-day calculation. Those skilled in the 

30 art will readily recognize that other, less accurate clocks such as GPS clocks can be 
used but the effective location resolution may be affected. 
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Differential Time of Arrival Svstem 



Figure 6 shows the basic principals of the Differential Time of Arriv al 
(DTOA) System. The example system is configured by at least lour receiving 
stations which operate at exactly the same clock. The transmitted signal is received 
5 by the receiving stations at T,. T 2 . T.. and T 4 respectively. The base stations arc 
located at (X Y ,). (X Y 3 ). (X Y ; ) and (X 4 , Y 4 ). The lour time 
measurements from the four base stations are resolved in non ambiguous fashion as 
follows: 

10 C - Velocity of Light. (1) 

X r Y, are the coordinates of the transmitter, and are resolved directly from 
those 3 equations- T, is a dummy variable representing the time of transmission. 



The four equations are dependent (rank 3). Therefore three linear equations could 

20 be from them (by subtracting row 1 from row 2. row 2 from row 3 and row 3 
from row 4): 

(Define: A P , , = P , - P, : where P = Oany parameter and CAT , , = A R , , ) (3) 

2AX I3 X, + 2AY I: Y t - 2CAR I: T, = X , : + Yf - X - Y : : - C : (T, : -T : : ) (3a) 

25 2AX „ X , + 2AY :3 Y, - 2CAR : , T, = X : 2 + Y : : - X r - Yr - C 2 (T : 2 -T, 2 ) (3b) 

2AX M X , + 2AY >4 Y t - 2CAR„ T, = X f + Y, 2 . X i - Y/ - C : (T\ 2 -T 4 2 ) (3c) 

Those three equations could be easily further reduced into two linear equations for 
X , , Y , dependent only upon AR ij. 
30 By dividing equation (3a) by A R,, and equation (3b) by A R : , and 

subtract the results and do the same for equations (3b) and (3c) we get the 
following equations. 



15 



(X,-X,) 2 + (Y , - YJ 2 
(X ; -X.) : + ( Y :; Y,) : 
(X 3 -X t ) 2 + ( Y 5 - Y,) 2 
(X 4 -X,) : + ( Y r Y,) : 



C 2 * ( T, - T, ) : 
C 2 *(T : -T ( ) : 
C 2 * (T,-T t ) : 
C 2 * (T^T,) 2 



(2) 
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(Define: I) M 2 = X , : + Y; - X - - Y, : ) (4) 
X, !A X i? /A R i: - A X r , /A R 2y \ + Y, <AY |: /A R i: - AY :; /AR : ,! = 

|(D I: 2 /AR i: ) - (D : f / AR r J - AR,,!/2 
X, |A X :? /A R : , - A X y4 /A R u | + Y, |AY\, /A R :: , - AY ;4 /A R ;4 [ = 
5 |(D : , : /A R :? ) - ( D : J I A R 3J ) - A R 24 ]/2 

The location could be resolved with only three bases, but it would 
result in an ambiguous solution which could be resolved by an a priori 
knowledge of the transmitter location. 

10 

Dual Frequency Differential Ranging Principle 
The differential ranging system is based upon transmitting two 
simultaneous frequencies for each hop, with an alternating difference that solves the 
ambiguity. The principal is described in Figure 8. The mathematical derivation of 
15 the differential ranging is based upon the phase difference that is created when two 
frequencies with a small difference in their wavelength are propagating over 
distance between the transmitter and the receiver. The two signals are represented 
as: 



20 A, (t) = A 0 * Cos (wt + 40 (5) 

A : (t) = A„ * Cos ((w-Aw)t + <|> 0 ) (6) 

The phase difference that is propagating through the distance R, 
between the transmitter and the first receiver is: 

25 

Acp, = wR, IC + * - (w - Aw) R, /C - 4> t( (7) 

The same applies also for the second receiver which is at a distance of R : 
from the transmitter: 

30 

A <p : = wR : /C + 4> - (w - Aw)R : /C - <j> „ (8) 
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Solving the iwo equations: 



AR i: = R, - R ; = (Atp, - A(p> )* CVAw (^) 



5 A <p, and A cp ? are measured by the digital receiver at each base station. Aw 

is known a priori. The basic algorithm is as follows: 

1. Compute the phase difference of the two frequencies at the first base, 
station. 

2. Compute the phase difference of the two frequencies at the second 
10 base station. 

3. Compute the differential range using formula (9) above. 



Ambiuuitv Resolution 
There are two types of ambiguities in the differential ranging system. 

15 First, there is ambiguity that is created if the signal is received by only three 
stations. In that case the hyperbolas are intersecting in two points. 
Mathematically, that ambiguity cannot be resolved unless we have an a priori 
information (for instance - one of the points is in the sea), tracking data (the true 
intersection will move on roads) or other means such as Interferometer Direction 

20 Finding. The second ambiguity is the ambiguity in estimating AR. The reason for 
the ambiguity occurrence in AR is apparent from formula (9) which should be 
expressed as: 

AR 12 = R, - R : = (Acp, - Aq> : + 27in)* C/Aw (10) 

25 

where n represents the fact that the difference in wavelength between the two 
frequencies could be larger than I. Formulating it (10) at the first base station: 



Define 4 J fj = (0 - 2*n) (11) 
30 R, = ( n„ + V„ /2n) * X = ( n !3 + T i: !2n) *(>. + AM t = R, / C 
A„(t) = A 0 * Cos <w( n M + 4' n /2n) * )JC + + ; 
w*A / C = (w - Aw)*(X + A?.) = 2n 
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A,,(t) = A n * Cos J 2k n M + M'„ + $\ 

A i: (t) = A„ * Cos «\v - Aw) t + f 0 ) = A„ * Cos j 2n n i: + M' i: + 4> n | 
AO-, - 2tt (n n - n i; ) + T n - + + - f, = 2n (n, ,-n,0 + AM', + <|> - 

5 For the second base station: 

R : = ( n 2I + M' 21 /2n) * l = ( n :: + 4\ 2 /2n) + AX ); t = R : / C (12) 
A : ,(t) = A„ * Cos { 2nn :i + H' 2I + + j 

A 2: (t) = A u * Cos ((w-Aw)t + <() 0 ) = A, * Cos [2nn :: + + fc, J 
10 Acp : = 27i(n :r n 2: ) + S' 2I -M' :: + <J>-^ (l =27t(n 3r n 22 ) + AH' : + (M<. 

By subtracting the two measured phase difference we deduce: 

A <p, - A <p 2 = 2n (n n - n, 2 - n 2 , + n 2: ) + AM J , - AH' 2 = A R I: *Aw/C (13) 

15 

In order for AR, 2 to be non-ambiguous n n - n, : - n :i + n 2: has to be 
0. AR I2 < D. In Figure 7, the length of one side of a triangle is longer than the 
difference between the other two sides. Summarizing it: 

20 Assume Aq>, - Acp : < n (14) 

2*Af * AR, 2 / C < 7i 
Af <C/(2* AR 12 ) 
Af < C /(2* D) 

25 For D practically limited by 50 Km, Af has to be smaller than 3 

KHz. We chose to use resolution of 2.5 KHz. The problem is that the accuracy of 
AR i: is inversely dependent upon Aw - the larger Aw is. the smaller the error in 
AR t2 . 



30 Propagated Wave Form 

The propagated wave form has to resolve two contradicting issues. 
It has a resolution of 2.5 KHz and it has a maximal difference in frequency in 
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order to reduce ranging error. The proposed wave form suggests that the frequency 
difference Aw will be an integer multiplication of the non ambiguous frequency 
<Sw: 

5 Aw, = i * k 'I c i c N|. N = Number of hops. (15) 

The phase difference is incrementing theoretically from hop to hop by ficp. not 
necessarily at the same order. 

10 6cps8w*AR/C (16) 

Acp, = Aw> * AR / C - i * 5w * AR / C = i*5cp 

The evaluation of AR is done by computing fkp which has no 
ambiguities. The proposed detection algorithm is as follows: 
15 1 . Create a complex vector of V, = { V (l Cos (Acp , V 0 Sin (Acp , )! . 

Those values are the quadrate output of the Acp , calculation and Acp , 
does not have to be calculated explicitly. 

2. Perform FFT ( V k ). 

3. Find the peak of the FFT and compute its exact "frequency." The 
20 "frequency" is a measure of the phase increments, and provides 

excellent estimation of 5'p. 

The algorithm automatically improves the signal to noise ratio by the 
factor of the number of hops. There is no dependency between the hop frequency 
and the frequency difference of the number of hops. The process is actually 

25 sampling in the frequency domain. The form of transformation from the frequency 
domain to the time domain is IFFT (Inverse Fast Fourier Transform) rather than 
FFT (Fast Fourier Transform). It has no effect on the results other than 
interpretation. The starting frequency does not have to be fiw or 0 but rather any 
frequency of our choice. It is equivalent to sampling in the time domain which 

30 does not start at T 0 = 0 which has no effect on the result of the FFT other than a 
phase shift. 
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Adaptation of Standard NexNct™ to Differential Ranuinu 
The overall system throughput could be reduced due to two impacts, 
hirst, the power of each bit is reduced by at least 2. due to the fact that the power 
is divided between 2 frequencies. In fact the power will be further reduced due to 
5 the non-linearity of the power amplifier which creates spurious frequencies that will 
consume part of the power. Second, the transmission of two frequencies requires a 
new frequency plan. The frequency resolution is reduced to 2.5 KHz (from 7.5 
KHz) and a special care must be taken to align the hopping with the transmission 
of the two frequencies. 
10 The creation of the two frequencies is done using a regular mixer 

which will multiply the local oscillator by a sine wave in Af/2. The result will be: 

Cos (2nf LO t + <t>, 0 ) * Cos ( 2n{Af/2}t + * 0 + <J> m(Kj ) = (17) 
0.5*Cos (2*{f UJ -+ Af/2}t + {<|> m + +„} + <|> moJ ) 
15 + 0.5*Cos(27i{f, o - Af/2}t + «> m - +„} - <t> mod ) 

fc^ s£ ±71/2 (BPSK) 

As is seen, the power is distributed between the two signals. The 
phase modulation is applied through the differential modulating frequency. Both 
20 frequencies are phase modulated in inverse which has no significant since the 
demodulation is differential. The frequency difference is Af. 

Design of the Remote Mobile Unit 
In Figure 9. the modulating frequency is generated by the 

25 microprocessor 801 from a port called "TCMP" (time compare) which is capable of 
generating a square wave pulse train. The pulse train is modulated by 0 : or 180' 
using a digital XOR gate 802. The modulated square wave is passed through a 
low-pass (or band-pass) filter 803 which clips its harmonics and produces a sine 
wave which produces the two frequencies in the mixer. The ranging transmitter of 

30 Figure 9 is designed to avoid carrier rejection and high order harmonics. The 
power amplifier performance is selected to be extremely linear. This results in 
avoidance of leakage of energy to other channels. 
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The ll : circuit for the ranging system is shown in Figure 10. The 
local oscillator downconvcrts the IF signals from 70 ± 0.65 Ml Iz to 21.4 ± 0.01 S 
MHz in 36 KHz channels. Two separate converters are used to operate in parallel. 
In order to ensure phase continuity between hops (the phase difference between the 
5 two frequencies must remain between hops) we must base the double loop 
synthesizer on DDS. (Direct Digital Synthesizer). 

The process of finding the phase difference between the two 
frequencies is based on convolution of the two signals: 

10 A<p = arg( S, ® S : + ) (18) 

where arg(*) is defined as the argument of the complex variable (•) The process is 
identical to the process of interferometer DF. The convolution is done on a 
complete hop (45 msec) after the data extraction and the necessary phase reversal. 

15 The phase data of each hop is transferred to the central station. The process of 
differential phase comparison is done in the central station. The location 
processing includes differential ranging computation as well as a combination of 
differential ranging and direction finding. In operation, differential ranging 
requires that periodic calibration be performed in real time. The accuracy of the 

20 system will be in the order of 300 nanoseconds, which mean that timing error in 
the range of 50 nanoseconds must be compensated. The most cost effective way to 
do it is to set up several transmitters in predefined places and to calibrate the 
timing offset accordingly. Assuming a short term drift of 10 ,( ' (high performance 
oscillator), the time between calibrations will be: 

25 T = 50* 10' 9 / 10* 1,1 = 500 seconds ( l ( » 

The present system supports remote activation of the calibration 
transmitters. 

30 Accuracy of DTQA 

The accuracy of the DTOA system is determined by comparing it to 
a regular triangulation direction measurement system. The geometric shape 
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representing the points that have the properly of the same difference lime of arrival 
relative to a couple of bases is a hyperbola as shown in Figure 7. According to 
trigonometric law of cosines: 

5 AR = R : -R, = {R : + (D/2)) : + RD Cos (G)) ,: -{R : +(D/2) : - R D Cos(G) J ,: 
(2R : + 2(D/2) 2 - AR 2 ) 2 = 4 [ ( R : + (D/2) 2 ) : -((RD Cos (G) } 3 ] - 

8R 2 (D/2) : - 4R 2 AR 2 - 4(d/2) : AR : = 8R : fD/2) : - 4 j ( R D Cos (G) \ : 
Cos (G) = (AR/D) ( 1 + D 2 /4R 2 V 2 for (AR/R) « I. (20) 

10 In order to simplify the evaluation of the inaccuracy of the location 

due to the inaccuracy of AR estimation, the inaccuracy of the asymptote is 
calculated. This is valid when D/2R«1. Assuming D/2R«1, the asymptote 
formula is: 

15 Cos (9 A ) = AR/D (21) 

The asymptote has an ambiguity of ±9 due to its being an even function. 
Estimating the error in G due to the error in AR is: 

20 a(0 A ) = a( AR) / D * Sin(0 A ) (22) 

For an example, let E(AR) - 50m. D=20 Km, E(G A ) is 0.14'. It 
could be shown that the hyperbola error is very close to the asymptote error. The 
error at the transmitter location could be evaluated as: 

25 

a(R) = Ra(9) = o( AR) / {(D/R) Sin(9)} (23) 

Note that the error in the estimation of the DTOA is multiplied by a 
geometric factor (1/ (D/R) Sin(9) ). The constellation of the two receiving bases 
30 and the transmitter could be scaled up or down and the effect on the errors is the 
same. The term (1/ (D/R) Sin (9) ) is recognized as GDOP. 
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GDOP = (D/IO Sin(U) (24) 
where CiDOP = geometric delusion of properties. 

Hrror Sources in DilTerenlial Ranninu 
5 External Thermal Noise is one source of error in the present 

invention. l ; rom equation (9) follows equation (25): 

a(AR i: ) = a( A<p i: )* C/Aw (25) 

10 For signal to noise ratio (SNR) above 10 dB the system uses the 

approximation (see Figure 11): 

a(Acp) = 1/ (2 SNR) r: : ..a ( A<p i: ) = 1/ (SNR) : : (26) 

15 The error in Atp is dependent upon the algorithm. In order to simplify the 

process, in the alternative, the following is an algorithm for the extraction of the 
frequency out of the FFT. Find the peak of the FFT and its strongest neighbor. 
The value of these peaks are marked as P and P 1 . The following is the amplitude 
distribution is a sine wave in FFT with boxcar window (signal power equals A„ 

20 which is a function of the signal to noise ratio induced by Acp ). 

A - A () * Sin (it - 6q) (27) 
N * Acp - 2rck + 6q 



2} 



30 



Since Sin (6c) = Sin ( n - 5q) the system deduces: 

A7A = 8q/(7t-5c) (28) 
8c = nA7(A + A') 

In order to evaluate the accuracy of the estimation, a point is chosen 
which is the worst regarding the power distribution in the cells of the FFT. At this 
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point A = A". Assume that ihe noise is white i.e. aA = a/Y - vn. The error in 6c 

is: 

a(5c) = 7i |aA( A : + A ,2 ) ,: /(A + A') : ! = \n/2\2\ * jaA/AJ (2<)) 
5 A « A „ * Sin ( n/2 ) = A u *2/n 

SNR,, = (A ( / aA) 2 
o( 5c ) = {tt/2} : /(2 * SNR 0 ) 12 

Since the input noise is distributed over N cells of the FFT: 

SNR„ = SNR j * N (30) 

SNRj is the input signal to noise ratio. Finally: 
a(AR l2 ) * a( Aq> 12 ) * C/Aw = ( n/8 ) * C / (5f*N 3/3 {2 * SNR, ] ,;2 ) (31) 
1 5 where C = 300,000,000 m/s 

N = 53 

Sf = 2500 Hz, 

SNRj = 10 dB = 10 (integrating over a complete hop) 
20 And the final result of the estimation of AR, : is: 

a(AR 13 ) * 27 m (32) 

Multinath as a Source of Error 

25 The effect of multipath is very different from the effect of the noise 

in two ways. First, the multipath is not necessarily independent, thus the averaging 
that is done through the FFT will not provide an improvement of N 3 : . Second, the 
signal to reflection ratio is rarely in the order of 10 dB. Several simulations were 
done in order to evaluate the effect of the multipath and the results showed that the 

30 frequency hopping scheme randomizes the effect of the multipath. 
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Dopplcr Shift as a Source of Error 
The Dopplcr effect is only on ihe frequency difference. In 1 MHz 
difference we have: 

Afdmav = V (mv /A = V m:iv Af/C 

5 Where: 

Af d(im = Maximum Dopplcr shift 

V mav = Maximum radial velocity (assume 40 m/s) 

Af = Maximum frequency difference (up to 1 MHz) 

10 Result in: 

Af dmax = 0.13Hz (34) 

Which is negligible since it causes a bias of 0.006 radians in each 
15 hop. which results in: 

EfAR™} = E{ A<p }* C/A\v = E{(Ac/N)i * C/Aw = 13.6 m 
E(p) - Error in parameter estimation 

20 RMU Frequency Inaccuracy as a Internal Source of Error 

The modulation of the dual frequency in the Remote Mobile Units 
(RMU) is done with a low cost oscillator which has a possible drift of 100 PPM. 
The error causes a direct error of 0.0001 in the differential ranging (sec 32) which 
is an error of several meters in the worst case. Timing is another source of internal 

25 error. The most difficult bias error to overcome is the liming error between 

receivers, h is obvious that even an atomic reference clock with an accuracy of 10" 
L results in 50 nanosecond error within 1.5 hours. The system design is based on 
calibration transmitters in predetermined places. 

Internal Phase and frequency Errors is another source of internal 

30 error. Differential frequency errors could be practically eliminated by using a 
TCXO (temperature compensated crystal oscillator) with I PPM accuracy. Phase 
errors could not easily be eliminated internally. But. if the frequency generators 
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arc phase continuous, the phase error creates a differentia! ranging bias error (see 
equation T32) which has to calibrate continuously by the calibration transmitters. 
Note (he possible group delay at the II-' crystal filters. 

Wave Form 

The transmitted wave form is a dual frequency transmission with a 
d.fference of minimum of 200 KHz to 350 KHz. with resolution of 2.5 KHz. 

The Two-Fre quencv Differential Rnno. Ai„„ ri , hm 
1 0 The following description is the complete mathematical analysis of 

the two frequency differential ranging algorithm taking into account all known 
sources of error. Let l denote the true time for the /th hop. where t = 0 
corresponds to the beginning of the /th hop. In the following equations, a dot 
above the variable indicates its derivative. The transmitter frequency at the 
1 5 beginning of the /th hop is given by 



L,=U, +«/, + 6f,(i- \)T. (T)) 

where 

./„., is the basic frequency of the /"th hop; 



20 



25 



where 
30 • 



6/ x is the transmitter frequency bias: 
5/ is the transmitter frequency drift; 
T is the time interval between hops. 

The two instantaneous phases of the transmitter frequencies are 

= (f,, + L, U + 0.56 f v ? + Kn + ^ + /J( (/) (T2aj 
" (L -L, )t + 0.56 f v v + . Kii + , h (t) (T2h) 



L, is the frequency of the DSB modulation at the /th hop; 
4>,„ is the initial carrier phase at the /th hop; 

is the initial phase of the DSB modulation at the /th hop: 
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», ID is Ihc phase noise of the transmitter oscillator. 
The instantaneous ranges at the /th hop are 

ID = '\„ + v u , \i + ( i . \yr\ + 0.5a,,, \i + (i - |;71 ? 

5 = k« + >v» r/- ur+ 0.5^ „ r/ - \fr : \ 

+ K« + <w- urj/ + o.5«„/ : 

= + »V, / + 0.5a lf) r. 

where 

A- = 1 .2 is the receiver number; 
10 • r u , is the range between the transmitter and the jfcth receiver at the 

beginning of the first hop; 

v K * is the velocity along LOS between the transmitter and the *th receiver 
at the beginning of the first hop; 

a kM is the acceleration along LOS between the transmitter and the *th 
1 5 receiver at the beginning of the first hop: 

i\, is the range between the transmitter and the *th receiver at the 
beginning of the /'th hop; 

»V, is the velocity along LOS between the transmitter and the Jtth receiver 
at the beginning of the /th hop. 
20 The two received signals at the kih receiver arc 

h , (i) = 2 ' - exp{/27i<t,,„ (i . * 0 /yW }'«<,. (,j + „- (< U) (T4u) 
M (0 = 2-' - exp{/2n<|>„ ft . K,^(i))\m u , (tj - „■„. <„ . (T4h] 

25 where 

k„ is the inverse of the speed of light: 
»/,„(/; and m tfl (ij are the multipaths of the two frequencies: 
»••„, (i) and wjt) are additive white noises for the two frequencies. 

30 ,n P ract 'ce. the two frequencies are not received separately, but their 

sum <„ (0 + u (I) is received. However, since the two frequencies are separated 
by 7.5 KHz at least, which is much above the bandwidth of each frequency, we 
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assume lhat ihcy can be separated by preiiltering. This also explains the factor 2' 
'\ namely, lo rcllcel the fael that the total transmitter power is divided between the 
two frequencies. 

The Alh receiver generates two reference signals whose instantaneous 

5 phases arc 

(0 = (U+LJV -0 + 6 t rk ) + 4> lM (T5</) 
(0 = (/'„,- LJO '0 + 5 in) + - .(T5b) 

10 where 

9 is the delay, with respect to the beginning of the /th hop. at which 
demodulation starts: 

6 i tX is the time bias of the kih receiver. 

<|> lU , is the initial phase of the reference signal of kih receiver at the /th hop. 
15 Note that the receivers are assumed to have "ideal" frequencies (zero 

frequency deviations), because the reference signals arc synthesized digitally, so /,,, 
and /„, , are just numbers in the computer. 

The demodulated signals at the receivers are 

20 d ktl (V = cxp {j2n^,„ (t)\ \ n (0 (T6a) 

d kh (() = exp (/2tt<|u (0) \ h (0 . (T6h) 

The equations derived so far are used in the simulation. In order to 
derive the estimation let us now ignore the noise. Also express the multipalh in 
25 terms of its magnitude and phase, that is. 

»h t > (0 = A to 0) exp U2n\ Vtu (i)} (11a) 

m kh (0 = A ih (0 exp {/2n\|/ M (t)/. (T7h) 

Then, 

30 d ka (t) = 2- V2 A ku (U exp U2n (!) - - *„/•, flj) - ! (T8a) 

tt kh (l) = 2*'- w exp {/2tc li<„ (U - <*>,,, 0 - *V, ^ - (T8/» 
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The instantaneous phase of each demodulated siunal contains a 
quadratic term, a linear term, a constant term, and a random term. Denote the 
coefficients of the linear terms by / ft „ and f M respectively, h is easy to verify that 
the linear terms are approximately given by 

5 

./;«„= «„v„(/;, +/,,)- bj\- ft f v ti - 1)7- (T9«) 
U = Ko v 0 C/i, - ./„,, ) - &/,- S i; (i - 1 )T. (T9b) 

Suppose we estimate the linear terms (using a Discrete Fourier 
10 Transform with interpolation), and use the estimates to form 

«, =/d,„ +/,/,* - In,, 'Inn * K oL., ( v i., ' »V, ) * K„f n p\ , - v,, ). (T10) 

The use of a , is explained below. 
1 5 Next let us form the signal 

iL (0 = dt.atfu, (0 = 0.5A U (t) A Vh (t) 
exp{/'2Tr {<t> jh (f - k„ r t (l)) - <f> M (i - k„ r t (I)) + 4> tlM (0 - <*> rlh (0 + 

K (Tin 

20 

We have 

*«* 0 - k„ r t (0) - <t> M (i - k„ r t (()) = -2f m , (i - k„ i\ (i» - 2<f>,„., 

= 1.L, k b '\, - 2f mi ( 1 - k„ v L , )i + f nu K 0 « M r - 2<*>,„, (T 1 2 ) 

25 4u, W - W = V m , ' - U,,, 0 + 2f m , 6/ r< . (T131 

So. 



& W = 0.5/1,,, (t)A kh (l) 
30 exp{/2ir [2/„,yc,/Y +2/„,yc,^^ H"14) 
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The maximum value of 2f mt K„vj is about .*H()"\ and the maximum value of /,'„,*„ 
a kU r is about 2*1 0 \ so these two terms can be neglected, leaving 

* 0.5/^,//M„//ycxp{/27r I2/ m ,KVo - 2^ r 2/ w ./»+2/^^ r ^^'^-^|l ■ fTI5) 

Assuming thai W and iw w M are strongly correlated, we have ^ Xw * $ L/t , Let 
0 < n < A' - 1 } be the sampling instants of the /th hop and form 

(T16) 



p lk ^(2n)' u angle 



N-l 

E i k U B ) 

73 = 0 



= 2 f mt s K c r k4 d -2<^ , -2 4, ,6 + 2 f ^ ,6 t , , 



10 and 



ft = ft, - ft : - 2f„ u K tl (r u - + 24,(41,, - «/ r2 ). (T17) 

The bandwidth of the signals djt) and dj\) is limited to about 200 
15 Hz (twice the maximum Doppler frequency). Therefore they should be low-pass 
filtered to this value and decimated to a sampling frequency of about 200 Hz. 
With the present value of T (50 milliseconds), we get 10 samples per hop. If the 
bit rate is also 10 per second, the range estimation algorithm can be conveniently 
combined with the bit-detection algorithm, since the phase of each bit is needed for 
20 both purposes. The simulation does not include the filtering, but the standard 
deviation of the additive noise is computed in agreement with the assumed 
bandwidth. 

Finally use the estimates a, from equation T10 to compute 
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i-l 

if = w - vJL a ' 

f=0 /,,f 

i-\ 

- 2f m fi 9 (r,, - r : ,> - 2/„,,k« £ (V " " ^ } T + 2J » JK ' ' 8/ ' = 1 

* 24,k 0 (/ m , - r Wl ) + 2^,(8/,, -8/ r3 ). (T18 ) 

10 As we see, the role of a, is to compensate for the motion during the 

hops and to estimate the initial range difference. Motion compensation depends on 
the receivers carrier frequencies being identical. If they cannot be made identical, 
this compensation will not be performed, and 0, will be used for subsequent 
processing instead of y t . 

] 5 The final stage of the algorithm is the estimation of the range 

difference from {y # , 0 < / < 1 - 1 J, (or from { P,, 0 < / < I - 1 } . as discussed 
below). Let £(i) be the ordering of the hops such that /,, , (i) = /„„„ + t f tnU .0 < f 
< / -I. Let r| f = exp' t /27tY,,„i. Then, except for the receiver synchronization 
errors 6 t rk . r| f is given by 

20 

x\ ( = exp{/47t Ko (r I0 - r 30 ) {f ttun + 2 f m(i )). (T19) 

Therefore. 2f„ fU K fl (r Ul - r 2ja ) can be estimated by performing a Discrete Fourier 
Transform with interpolation on the sequence {r| 0 < f < / - 1 J . 

25 If we use /?, instead of 7,, the sequence of r\ t will not be a pure 

sinusoid (even in the ideal case), due to the transmitter motion. However, if 
is sufficiently random, the phase deviation seen by (using a Discrete Fourier 
Transform with interpolation) will be random, so this should not have much effect 
on the estimate. For example, if the speed is 30 m/s. the transmitter motion during 

30 the entire sequence of hops is no more than 120 m. With /„„„ + f w U (/ - I ) = 60 
KHz. this translates to about 0.05 # 2n radians phase error. It is therefore expected 
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lhat ihc Discreic I'ouricr Transform wilh interpolation procedure will give a good 
estimate of the average range, or the range at the (//2)lh hop. 

The Multipath Model 

5 Multipath is basically a linear phenomenon, that is. the received 

signal is a superposition of several signals, corresponding to several propagation 
paths from the transmitter to the receiver. One of these is the direct path, whose 
impulse response is assumed to be an impulse, with a gain factor a and a zero 
delay (the nominal propagation delay to the receiver is immaterial, and can be 

10 ignored). All the rest are reflected paths. Each reflected path is characterized by a 
random gain and a random delay. These random parameters vary over time in 
general, so the total response at time / to an impulse transmitted at time / 0 is given 
by 

h(t, t 0 ) =ad(t-t 0 ) +Lcc n (t 0 )6(t-t 0 -T n (t 0 )) . [T20) 

15 The case a = 0 is called Rayleigh fading, while the case a > 0 is called Rice fading. 
In the Rayleigh case there is no direct propagation, only reflections. This happens 
when there is no line of sight between the transmitter and the receiver, and the 
wavelength is too short for diffraction to have an effect. 

When the transmitted signal is sOJ. the received signal is given by 

20 the convolution of sfl) with h(u ( u \ that is. 

r (t) = as(t) +Ea r (t)s(t-x n (t)) . (T21) 

If we now assume the existence of infinitely many reflected paths, each having an 
infinitesimal contribution. (T2I) becomes 

r (t) =as (t) + f°°a(x t t)s(t-T)dr, (T22) 
Jo 

25 where ol(tX) is the infinitesimal contribution at time / of the reflection delayed by 
t. The statistical properties of afr,t) will be discussed later. 
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Now Id sit) be a C\V signal, namely. s(n - cxp( [ilnfy] . Then 
r (c) ~as(ti * y~a (x, t)e>ipii2nf 0 {c-r) ) cx -cxp \ i2n f c C) \a*x(t)\ , i T3J ) 

where 

m(t) - Fa (x, C)exp{-j2nf 0 z) dr. . ( T24 ) 

Jo 

As we see from (T24). the effect of multipath on a CW signal is 
multiplicative, and it depends on the function of a(r.t). Here we are interested in 
5 short-term multipath in urban areas and mobile transmitters (or receivers). A 
common model for a(rj) in such scenarios is as a Gaussian process, stationary and 
narrowband in /. and nonstationary but uncorrelated in r. Thus we assume 
Eol(tJ) = 0, 

r(x, , x 2 : /, , ( 2 ) = Ea(r } ,/,)a(x : • '2) = r , Ci - '2) r , ( t i) 5 ( t i ' T :)- (T25) 
10 The temporal correlation T, (7, - t 2 ) depends mainly on the physical 

motion of the transmitter. In particular, if the doppler frequency of the transmitter 
is f thppfi . r (corresponding to the wavelength and the transmitter speed), then the 
spectral density of T, (t x - / 2 ) is limited to ±f tUipvk . r . The shape of the spectral 
density is difficult to estimate, since it depends on the given scenario, mainly on 
15 the spatial distribution of the reflections. In the simulation we used a flat spectrum 
in the range ±f thppfcr for simplicity and "universality." 

The function Y x (x,) describes the relative multipath power at 
different delays. It is also scenario-dependent and difficult to estimate. However, 
based on the graphs known in the art, we assume a functional dependence of the 
20 form 

f f (T^exp{-T,/T mp ), (T261 
where t is a scenario-dependent parameter. Note that when I\ r r, ) is plotted in 
dB. it shows a decreasing straight line. 

It is interesting to compute the multipath correlation for both 
25 different times and different frequencies. Let m{t x . /,) be the multipath at 
frequency./, and time f, . and similarly for m(t z . / : ). Then 
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r.,U., £,;/,,/,) -Em(t,. £.)m*(t :l f .) 

= j^j^r, ( t , - c 2 ) T, ( t , ) 6 ( t ! - 1 2 ) exp { - j'2 k ( f ,-c , - f 2 t 2 ) } dT , dt 2 



10 0 



= fy t ( tj - 1 2 ) TJtj) exp { - j - 2n ( f j -f ,) t, ) } eft, 



l+j'2n ifj-f 2 ) x 



mp 
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Table 1: Simulation Results lor the Two I'raiucncv Model 
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Table 1 shows the results of simulation of the present system with 
the given parameters. The frequency f^ K is selected as described above as/„ wl + 
Lo (/ - 1) SNR is the S| gn al to noise ratio, MP = 0 indicates Rayleigh fading 
15 and MP = 1 indicates Rice Fading. STD is the standard deviation (one a), and the 
resolving power of the algorithm is given as minimums (Min) and maximums 
(Max). 

20 Although specific embodiments have been illustrated and described herein, 

it will be appreciated by those of ordinary skill in the art that any arrangement 
which is calculated to achieve the same purpose may be substituted for the specific 
embodiment shown. This application is intended to cover any adaptations or 
variations of the present invention. 
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CLAIMS 

1. A mobile transmitter location system, comprising: 

a mobile transmitter having a unique identifier, and having a 
transmitter for simullanctuisly sending a first frequency-hopped spread spectrum 
5 carrier signal and a second frequency-hopped spread spectrum carrier signal, at 
least one of which containing the unique identifier, the first frequency-hopped 
spread spectrum carrier signal having a different frequency than the second 
frequency-hopped spread spectrum carrier signal; 

a first base station having an antenna operable for receiving the first 
10 and second frequency-hopped spread spectrum carrier signals and capable of 

determining a first time of arrival of the frequency-hopped spread spectrum carrier 
signals; 

a second base station also having an antenna operable for receiving 
the first and second frequency-hopped spread spectrum carrier signals and capable 
15 of determining a second time of arrival of the frequency-hopped spread spectrum 
carrier signals; and 

a central station connected to the first base station and the second 
base station, and operable for: 

1. ) receiving the first time of arrival of the first and second 
20 frequency-hopped spread spectrum carrier signals from the first base 

station; 

2. ) receiving the second time of arrival of the first and 
second frequency-hopped spread spectrum carrier signals from the 
second base station; 

25 3.) generating a first hyperbolic plot based on the differential 

time of arrival of the first and second frequency-hopped spread 
spectrum carrier signals at the first and second base stations; and 

4.) using a third information source to determine the location 
of the transmitter along the first hyperbolic plot. 

30 

2. The mobile location system according to claim 1 wherein a reference station 
operates periodically to transmit a reference first and second frequency-hopped 
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spread spectrum carrier signal from a known location to calibrate the first and 
second base stations and the central station. 

3. The mobile location system according to claim 1 wherein the first base 

5 station and the second base station both include atomic clocks to synchronize the 
lime of arrival calculations. 

4. The mobile location system according to claim 1 wherein the third 
information source is a geographic map of the area around the base stations which 

10 provides information on the probable location of the transmitter along the first 
hyperbolic plot. 

5. The mobile location system according to claim 1 wherein the third 
information source is a third base station which calculates the time of arrival of the 

15 first and second frequency-hopped spread spectrum carrier signals and wherein the 
base station generates a second hyperbolic plot to intersect the first hyperbolic plot 
to identify the location of the transmitter. 

6. The mobile location system according to claim 1 wherein the transmitter 
transmits a plurality of hops each having a different first and second frequency- 
hopped spread spectrum carrier signals and the time of arrival is determined by 

a. ) calculating the phase difference between the first and second frequency- 
hopped spread spectrum carrier signals for each of the plurality of hops to produce 
a plurality of phase differentials: and 

b. ) performing an Inverse Fourier Transform on the plurality of phase 
differentials to produces a time of arrival. 

7. A method of locating a mobile transmitter location, comprising the 
steps of: 

30 receiving at a first base station a dual radio carrier having two 

frequencies and calculating therefrom a first time of arrival of the dual carrier: 
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receiving at a second base station a dual radio 
carrier having two frequencies and calculating therefrom 
a second time of arrival of the dual carrier; 

comparing the first time of arrival with the second 
5 time of arrival; and 

determining the location of the origin of the dual 
radio carrier by hyperbolic intersection location. 

8. A mobile transmitter having a unique identifier for 
10 use in the method of claim 7 comprising a transmitter for 

simultaneously sending a first frequency-hopped spread 
spectrum carrier signal and a second frequency-hopped 
spread spectrum carrier signal, at least one of which 
containing the unique identifier, the first frequency- 
15 hopped spread spectrum carrier signal having a different 
frequency than the second frequency-hopped spread 
spectrum carrier signal. 

9. A base station for use in the method of claim 7 
20 comprising an antenna for receiving the dual radio 

carrier and means for determining the first or second 
time of arrival of the dual carrier. 

10. A mobile transmitter location system comprising: 
25 a first base station having an antenna for receiving 

a dual radio carrier having two frequencies and means to 
calculate therefrom a first time of arrival of the dual 



WO 96/25673 



PCT/GB96/00270 



46 

carrier; 

a second base station having an antenna for 
receiving a dual radio carrier having two frequencies and 
means to calculate therefrom a second time of arrival of 
5 the dual carrier; 

comparing means for comparing the first time of 
arrival with the second time of arrival; and 

location determination means for determining the 
location of the dual radio carrier by hyperbolic 
10 intersection location . 
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